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ABSTRACT 

We have obtained a Ha position-velocity cube from Fabry-Perot interfero- 
metric observations of the HH 110 flow. We analyze the results in terms of 
anisotropic wavelet transforms, from which we derive the spatial distribution of 
the knots as well as their characteristic sizes (along and across the outflow axis). 
We then study the spatial behaviour of the line width and the central radial 
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velocity. The results are interpreted in terms of a simple "mean flow+turbulent 
eddy" jet/wake model. We find that most of the observed kinematics appear 
to be a direct result of the mean flow, on which are superposed low amplitude 
(~ 35 km s~^) turbulent velocities. 

Subject headings: ISM: Herbig-Haro objects - ISM: jets and outflows - ISM: 
kinematics and dynamics - ISM: individual (HH 110) 

1. Introduction 

The HH 110 Herbig-Haro jet was discovered by Reipurth & Olberg (1991), and it 
presents a complex structure of aligned knots, which opens out with an unusually large 
angle of ~ 10°. While HH objects such as HH 34 or HH 111 show aligned knots with more 
or less organized, arc-like shapes, HH 110 has a more chaotic structure, more reminiscent of 
a turbulent flow. 

A possible explanation for the fact that HH 110 is different from other HH jets has 
been proposed by Reipurth, Raga & Heathcote (1996), who suggested that HH 110 is the 
result of a deflection of the adjacent HH 270 jet (which travels approximately in an E-W 
direction) through a collision with a dense molecular cloud core. This interpretation was 
supported by the fact that while the source of the HH 270 jet is detected in the radio free- 
free continuum, no such source is detected along the HH 110 axis (Rodriguez et al. 1998). 
Also, infrared images show H2 emission concentrated on the W side of the HH 110 jet (Davis, 
Mundt & Eisloffel 1994; Noriega-Crespo et al. 1996), which would correspond to the region 
in which the jet is in contact with the dense cloud which produces the deflection. Against 
the jet/cloud collision scenario are the results of Choi (2001), who do not find molecular 
(HCO''") emission around the jet/cloud impact point, as might be expected. 

There are a series of papers which describe analytic and numerical models meant specif- 
ically for HH 110 as a jet/cloud impact flow (Raga & Canto 1995; de Gouveia dal Pino 1999; 
Hurka, Schmid-Burgk & Hardee 1999). The latest models of this object are the ones of Raga 
et al. (2002), who have carried out the most detailed simulations up to now of this object, 
and computed Ha and molecular hydrogen emission maps which agree with the observations 
in a quahtatively successful way. These authors also suggest that the point at which the HH 
270 jet is currently impacting the cloud might have moved westwards of the HH 110 axis (i. 
e., deeper into the cloud), which would explain the lack of HCO"'" emission in the crossing 
region of the HH 270 and HH 110 axes. 

In a recent paper, Riera et al. (2003) have discussed the kinematics and excitation along 



-3- 



and across the HH 110 outflow, as derived from long-slit spectra obtained with 5 sht positions 
(one along and four across the outflow axis), finding very complex structures. In the present 
paper, we discuss new Fabry-Perot data which cover only the Ha line, but give a full spatial 
coverage of the well resolved HH 110 outflow. The full spatial coverage of the position- 
velocity cube obtained from our data allows us to carry out a somewhat novel analysis (at 
least, in the context of HH objects). As the observed spatial and kinematical structures 
are very chaotic, and appear to defy a direct description, we have chosen to describe these 
structures through a wavelet analysis. A wavelet transform of an image gives information 
similar to the one that could be obtained from a more standard Fourier analysis, but also 
preserves the spatial information. 

In this way one can obtain, for example, the characteristic sizes (along and across the 
outflow axis) of the knots as a function of position along the jet. This kind of information 
could also be retrieved by carrying out flts to the intensity proflles of the knots (e. g., with 
anisotropic Gaussian functions, as has been done in the past for other objects by Raga, 
Mundt & Ray 1991), but such an approach is impractical for an object as complex as HH 
110, in which knots of different shapes coexist across the sections of the outflow. 

We therefore obtain a position-velocity cube from our Fabry-Perot observations (de- 
scribed in §2), and obtain velocity channel maps from which we compute an Ha image, 
as well as radial velocity and line width maps (see §3). We then carry out an anisotropic 
wavelet analysis of the Ha image, and determine the spatial scales present in the map as a 
function of position along the outflow (§4). We then study the variations of the line width 
and the radial velocity of the emission as a function of spatial scale and of position along 
and across the jet (§5). Finally, we interpret this kinematic information in terms of a "mean 
flow-|-turbulent eddy" model of a turbulent jet (or wake), and use the model to deduce 
the mean flow parameters that are relevant for the HH 110 jet (§6). The results are then 
summarized in §7. 



2. Observations 

The observations were obtained at the European Southern Observatory 3.6m telescope 
in January 1997. The Fabry-Perot instrument CIGALE ^ was used. It is composed of a 
focal reducer (bringing the original f/8 focal ratio of the Gassegrain focus to f/2), a scanning 
Fabry-Perot and an Image Photon Gounting System (IPGS). Table 1 contains the journal of 



^The instrument CIGALE (for CInematique des GALaxieEs) is a visiting instrument belonging to the 
Laboratoire d'Astrophysique de Marseille. 
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the observations. We have used a Fabry-Perot with an interference order p — 793 (380 km 
s"-*^) and a Finesse F = 11. The IPCS, with a time samphng of 1/50 second and zero readout 
noise, makes it possible to scan the interferometer rapidly, avoiding sky transparency, air- 
mass and seeing variation problems during the exposures and thus has several advantages 
over a CCD for this application. The exposure time per elementary channel was 5 seconds. 
Several observations of HH 110 have been obtained in Ha during 2 nights. The two sets of 
observations in the Ha line have been done using 2 different interference filters in order to 
check if the shape of the interference filter affects the shape of the Ha profiles. A 18300 s 
exposure was obtained with a Aq = 6563 A (central wavelength) and AA = 9 A (FWHM of 
transmission function) filter, and a 3600 s exposure was obtained with a Aq = 6564 A and 
AA = 17 A filter. 

The reduction of the data cubes was performed using the CIGALE software. The data 
reduction procedure has been extensively described in Amram et al. (1992, 1993, 1994, 1995, 
1996, 1998). Wavelength calibration was obtained by scanning the narrow Ne 6599 A hue 
under the same conditions as the observations. Velocities measured relative to the systemic 
velocity are very accurate, with an error of a fraction of a channel width (< 3 km s~^) over 
the whole field. 

The signal measured along the scanning sequence was separated into two parts: (1) an 
almost constant level produced by the continuum light in a narrow passband around the line 
and (2) a varying part produced by the line. The continuum level was taken to be the mean 
of the three faintest channels, to avoid channel noise effects. The integrated line flux map 
was obtained by integrating the monochromatic profile in each pixel. The velocity sampling 
was 16 km s"-*^. Strong OH night sky lines passing through the filters were subtracted by 
determining the level of emission from extended regions away from the HH object (Laval et 
al. 1987). 

3. The Ha channel maps, radial velocities and line widths 

Figure 1 shows the velocity channel maps obtained from the Fabry-Perot data of HH 
110. We detect Ha emission in hehocentric radial velocities ranging from —100 to -|-80 km 
s~^, with peak intensities in the —38.3 and —22.5 km s~^ channel maps. 

From the line profiles at each position, we have computed the velocity of the emission 
peak (with a parabolic fit) and the FWHM (with linear fits to the sides of the line profiles) , 
as well as the line flux (by adding the fluxes of the individual channels). The results are 
shown in Figure 2. 
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The integrated Ha map obtained from the present data looks similar to previously 
published Ha images of this object, showing a complex structure of knots that broadens 
out to the South. The radial velocities have low values (~ — > —60 km s~^), with more 
negative values to the South (see Fig. 2). The radial velocity structure across the jet is 
complex, showing variabilities which do not appear to show systematic trends. These results 
are qualitatively consistent with the radial velocities of this flow which have been previously 
obtained from long-sht spectroscopy (Reipurth et al. 1996; Riera et al. 2003). The FWHM 
of the Ha line ranges from 50 to 100 km s~^, and also shows a complex spatial structure (see 
Fig. 2). 

Because of the complexity of the observed structures, we have carried out an analysis 
based on wavelet transforms, which allows us to obtain a general description of the flow 
without attempting a somewhat hopeless "knot by knot" description of the observations. 

4. Wavelet analysis of the Ha image 

4.1. The anisotropic wavelet transforms 

Wavelet transforms have been used in different astrophysical contexts. For example. Gill 
& Henriksen (1990) have used wavelet transforms to describe the complex structures observed 
in CO radial velocities of molecular clouds. Obtaining a description of the complex structure 
of the HH 110 jet belongs to the same class of problem, i.e., identifying the characteristics 
of complex spatial structures without loosing the positional information. Therefore, we have 
applied a wavelet transform analysis to our Fabry-Perot data. 

A 2D wavelet transform analysis of the Ha image of HH 110 can be used to obtain the 
positions of all of the "knots" , as well as their characteristic sizes. This kind of procedure is 
described, e. g., in the book of Holschneider (1995). 

Once we have obtained the positions and the characteristic sizes of the Ha structures of 
the HH 110 jet from the wavelet spectrum, we calculate the spatial averages of the deviations 
of the line center velocity and if the line widths (averaging over the characteristic sizes of 
the structures) in order to study the turbulence in the HH 110 jet (sec §5). This procedure 
is completely equivalent to the one followed by Gill & Henriksen (1990). 

An interesting difference between our case and the "turbulent molecular cloud" study 
of Gill & Henriksen (1990) is that a jet flow has two clearly deflned directions : along and 
across the flow axis. Because of this, it is reasonable to expect that the characteristics of 
the structures along and across the flow axis might differ. To allow for this effect, we flrst 
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rotate the Ha image so that the flow axis is parallel to the ordinate. On this rotated image, 
we then carry out a decomposition in a basis of anisotropic wavelets, which have different 
sizes along and across the outflow axis. 

In particular, we have used a basis of "Mexican hat" wavelets of the form 

where r = [{x/axY + and and ay are the scale lengths of the wavelet along 

the X- and y-axes, respectively. After trying several different possibilities, we have chosen a 
C — (a^+a^)~^/^ normalization. It is possible to use other normalizations (e. g., C — y^o^o^), 
which lead to slightly different results, as the power is rearranged over different regions of 
the wavelet spectrum. However, the positions of the different peaks of the wavelet spectrum 
(see below) are not strongly affected by the choice of normalization. 

The process that is carried out in order to compute the wavelet transform is to choose a 
range for and ay (which are taken to have integer values of pixels). For all of the possible 
{ax I ay) pairs we then compute the convolutions 

Ta.,aAx,y)=J j I{x\y')g{r')dx'dy', (2) 

where r' = {[{x' — x) / a^]^ + [{y' — y) / ay^Y^'^ , and/(a;, y) is the observed map as a function of 
the pixel position (x, y) . These convolutions are calculated with a standard FFT algorithm. 

With this process, we then produce the transform maps Ta^^ay{x,y), corresponding to 
versions of the observed map I{x,y) which have been smoothed with the different ga,,,ay 
wavelets. These maps can be used as follows. 

We first fix a position y along the observed jet. Keeping this value of y fixed, we find 
the finite set of values (xfc, a^,- ^, a^.^) corresponding to the position across the jet and the 
spatial scales (across and along the jet) at which the wavelet transform has local maxima in 
the a^ and ay dimensions as well as in either the x or in the y dimensions. In this way, for 
each y we detect the positions Xk of the knots (which correspond to intensity maxima either 
along or across the jet axis), as well as their characteristic sizes {ax,k and ay,jfc). In practice, 
at many positions y along the jet one finds a single maximum, but in some positions up to 
three maxima are found, corresponding to different structures observed across the HH 110 
jet. 
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4.2. The Ha "image' 



The results obtained with the process described in §4.1 are shown in Figure 3. This 
figure shows the Ha image (see §3), which we have rotated by 13° so that the outflow is 
more or less parallel to the j/-axis. We have then convolved this map with a set of wavelets 
with 1 < < 20 pixels (with 0".91 per pixel) and 1 < < 20 pixels. The values of the 
peaks Xk of the wavelet transform obtained for each position y along the jet are shown in 
the right hand panel of Figure 3. Also shown (as error bars) are the values of ax,k and ay^k 
corresponding to each peak, which give an estimate of the characteristic sizes (across and 
along the jet, respectively) of the observed structures. 

In order to quantify the observed broadening of the jet, for each position y we compute 
the weighted mean of the x-spatial scale (perpendicular to the outflow axis) : 



Figure 4 shows < Oj; > as a function of position y along the HH 110 flow. In this flgure, 

y = corresponds to the position of the peak of knot A, and y grows moving South along 
the HH 110 axis. From knot A to the middle of knot C (y = — > 20"), < > grows from 
~ 3" to ~ 7". For y = 20" — > 70", < > increases more or less monotonically from ~ 7" 



dispersion (with < Gx > values ranging from ~ 10" to ~ 13"). For y > 90", < Gx > shows a 
wide dispersion, with values ranging from ~ 1" up to ~ 15". 

If we draw a straight line in Figure 4, passing through the peak < a^: > of knots A and 
L, we see that it forms an upper envelope to all of the values of < > in this region (which 
includes all knots from A to L). From this upper envelope, we deduce that the emitting 
region (from y = 0" ^ 70") broadens with a half-opening angle of 10°. 

Finally, in Figure 5 we show the Qjy ^ dbS db function of ttx^k for different regions along the 
HH 110 jet. We plotted together the points with values of y in flve different intervals along 
the jet beam (knots A-H, I-N, 0-Q, R1-R2 and S, corresponding to the y-intervals given in 
Table 2). We flnd that for knots A-H, the Uy^k (i- e., the spatial scales along the outflow 
axis) range from 1 to 18", while the ax,k (i- e., the spatial scales across the outflow axis) 
have values of up to 15". 

For knots I-N, the ax,k values are either smaller than 6" or larger than 12" . Interestingly, 
the intermediate scales, which are detected in other regions of the jet, do not appear as local 
maxima in the wavelet spectrum, and are therefore not picked up by the analysis described 
in §4.1. 




(3) 



to ~ 15", For y = 70' 



1" 



90", < ttx > shows a more or less constant value with a small 
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While basically no correlation is seen between the Uy^k and ax,k values in knots A-N, 
at larger distances from the source, these values appear to be correlated. For knots 0-Q, 
the points more or less fall on a ay^k/cix,k = 1-3 line, even though there is a wide scatter for 
0'x,k < 6". For knots R1-R2, the points roughly fall on ay^k/0"x,k = 0.5 line, with the exception 
of few points which show large scales along and across the beam of the jet. Therefore, we 
see a transition from structures that are elongated along the outflow axis (for knots 0-Q) 
to transversally oriented structures (for knots R1-R2). There are, however, some elongated 
structures at knots 0-Q, and some large structures both across and along the beam of the jet 
at the location of knots R1-R2 (see Fig. 5). Finally, for the S-shaped structure (at y > 140", 
which we have labeled S) both axially elongated and transversal structures appear to be 
present. 



5. The spatial distributions of the radial velocities and the line widths 

In order to describe the spatial dependence of the kinematical properties of the HH 110 
flow, we compute two different moments of the line profiles: 

^,^Jiv-V^^ (5) 

for each spatial pixel of the position- velocity cube. In these equations, v is the radial velocity, 
and Iv the intensity (at a fixed position x, y) of the successive channel maps. We carry out 
the integrals over the radial velocity Vr by replacing them with sums over the velocity channel 
maps. Vc (see equation 4) is the barycenter of the line profile (which from now we will call 
the "line center" radial velocity), and W (equation 5) is a second order moment that reflects 
the width of the line profile. 

We find that the value of W is quite strongly affected by the H 11 region (which is 
particularly bright in the surroundings of the faint S region, see figure 2). Therefore, we 
have subtracted this emission by fitting a linear dependence between the fluxes (seen in each 
channel map) in two rectangular regions (with a 10" width) on each side of the HH 110 jet. 

With the values of Vc and W computed for all positions (a;, y) on the plane of the sky, 
we then compute the following spatial averages: 

L W\x',y')I(x',y')dx'dy' 
I{x',y')dx'dy' ' 
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J,^^^V.{x\y')I{x\y')dx'dy' 

Is [Vc{x', y')- < K > {x, y)f I{x', y')dx'dy' 

where I{x',y') is the H« flux obtained from co-adding all of the channel maps. These 
integrals are carried out over areas Sa^^ay which are ellipses with central positions {x,y) and 
major/minor axes and ay corresponding to all of the wavelets that have been identifled 
as peaks of the wavelet spectrum (see §4.1 and figure 3). < >^/^ ix,y) corresponds to 
the line width spatially averaged over the ellipse Sa^^ay with a weight I{x',y') (see equation 
6). We then compute the average of the fine center velocities < Vc > within the ellipse, as 
well as the standard deviation < A^;^ >^/'^ of these velocities (see equations 7 and 8). Such 
deviations of the line center velocity, averaged over sizes chosen from a wavelet spectrum, 
have been previously used by Gill & Henriksen (1990) to study turbulence in molecular 
clouds. 

We then obtain values for < >^/^, < K > and < Av'^ >^/^ at all of the central 
points of the wavelets shown in Figure 3 (corresponding to peaks in the wavelet spectrum). 
In Figure 6, we show the line centers, widths and standard deviations as a function of position 
y along the jet (all of the points at different positions x across the jet beam and with different 
and plotted) . 

< > as a function of position y shows quite a wide scatter, and a general trend of 
more negative velocities from knot A (< Vc >~ —5 km s^^) to knots P-Q (< Vc >~ —35 km 
s~^). If we measure these radial velocities with respect to the surrounding molecular cloud 
(with a heliocentric radial velocity of +23 km s^^, see Reipurth & Olberg 1991) and correct 
for an orientation angle of 35° with respect to the plane of the sky (Riera et al. 2003), we 
then obtain that the (de-projected) jet velocity grows from ~ 50 to ~ 100 km s~^ from 
knot A out to knots P-Q. This is in good agreement with the results obtained from long-sht 
spectra by Riera et al. (2003). 

At larger distances (y > 120"), < Vc > grows, reaching ^ a,t y — 145", and then 
remaining with slightly negative values up to the end of the detected emission. We believe 
that this behaviour of < K > for y > 120" could be a result of changes in the orientation 
angle of the outflow with respect to the plane of the sky (because in this region, substantial 
side-to-side excursions of the jet beam are observed on the plane of the sky, see Figure 3), 
and might therefore not reflect a real slowing down of the jet flow. 

The standard deviation < Av'^ >^/^ of the line center velocity (see Figure 6) appears 
not to have any dependence on position y along the jet beam (showing values mostly within 
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the 7 to 15 km s"^ range) out to knot M (y ~ 75"). Beyond knot N {y > 90"), < Av'^ >^l'^ 
shows a more or less monotonia growth, reaching a value of ~ 25 km s~^ at y = 175". 

The average hne widths < W"^ >^/^ do not show a systematic trend with position along 
the jet, mostly staying within a range from 35 to 45 km s~^. This width of course includes 
the instrumental profile, which has a value of ~ 32 km s~^. 

An interesting feature is that many of the points in the < K >, < Av^ >^/^ and 

< W"^ >^/^ vs. y plots (see figure 3) appear to fall on more or less monotonic curves. This 
is somewhat surprising given the complex structure of the HH 110 outflow. As might be ex- 
pected, the points which fall on the monotonic curves are the spatial averages corresponding 
to the wavelets with larger sizes (a^;, ay ^ 10", see figures 4 and 5), which naturally give a 
smoother dependence of the average quantities as a function of position along the jet. 

The fact that the averages over sizes of ~ 10" give smooth curves as a function of 
position indicates that the kinematical properties of the flow have structures only in sizes 

< 10", which are superimposed on a smooth dependence on position along the jet axis. 
To some extent, this result is a quantification of the fact that while the Ha intensity map 
of HH 110 has complex, high contrast structures with many different spatial scales, the 
radial velocity and line width maps of this object only show small scale structures which are 
superimposed on a generally smooth distribution as a function of position along the jet (see 
figure 2). 

In Figure 7, wc plot the line centers, widths and standard deviations as a function of 
position X across the jet (all of the points with different Qx and ay are plotted). We plotted 
together the points with values of y in five different intervals along the jet beam (given in 
Table 2). 

< Vc > shows basically flat structures as a function of x (i. e., across the jet beam) in 
the regions of knots A to N. For the three regions at larger distances from the source (0-Q, 
R1-R2 and S), either growing or decreasing < Vc > vs. x trends are seen. The average hne 
width < W"^ >^/^ shows basically no trend as a function of x. 

For knots A-N, the dispersion of the line center < Av"^ >^/^ shows a central region 
(with \x\ < 6") with values ranging from zero to 12 km s"^, and two outer peaks (with 
= 6 — > 10") with larger velocity dispersions of up to ~ 20 km s~^. The central, low 
velocity dispersion region therefore has a diameter of ~ 12" (8 x 10^^ cm at 460 pc), and is 
surrounded by a high velocity dispersion envelope with an outer diameter of ~ 16" (1.1 x 10^^ 
cm). 

Finally, in Figure 8 we show the values of < A^;^ >^/^ and < W"^ >^/^ as a function of 
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the size Qx (across the jet beam) of the eUipses over which the averages have been computed. 
Again, we divide the points into the five chosen regions along the outflow axis. 

Prom Figure 8, it is clear that both < Ai>^ >^/^ and < W'^ >^/^ show very clear trends 
as a function of the size scale a^. The < W'^ >^/^ values are approximately constant as a 
function of a^, and show a decreasing scatter for larger ax- The < Av"^ >^/^ values also show 
a decreasing scatter as a function of Qx, but this decrease in scatter is superimposed on a 
shallow trend of increasing < Av'^ >^/^ vs. a-^. 

Interestingly, in this representation (of < Ai>^ >^/^ and < W'^ >^/^ vs. a-^), the fine 
widths and central velocity variations show a quite systematic behaviour. Because of this, 
we choose this particular representation in order to try to interpret the kinematical charac- 
teristics of the HH 110 jet. This is done in the following section. 

6. A turbulent wake (or jet) model 

It is not inmediately clear what the observed line widths and radial velocity dispersions 
as a function of spatial scale are actually telling us. In order to interpret these results, we 
use a simple, analytic model (described in Canto, Raga & Riera 2003), which allows at least 
a partial understanding of the implications of the observations. This model is as follows. 

It is known from laboratory experiments that the mean flow velocity of a turbulent jet 
or wake has a centrally peaked cross section. Superimposed on this flow are the turbulent 
fluctuations, which cannot be modeled in detail. 

HH 110 has a complex structure, probably formed by irregularly shaped shocks which 
give rise to the observed emission fine knots. The observed spatially resolved line proflles 
are then the result of a mean motion along the jet flow (the "mean flow") combined with 
the scattered motions resulting from the non-planar shocks (the "turbulent velocities", in 
analogy to the standard description of a laboratory jet). The distribution of the radial 
velocity (projected along the line of sight) over the cross section of the jet will then be 
determined by the appropriately projected distributions of the mean and turbulent velocities. 

As a simple parametrization of the mean flow velocity, we assume that it is directed 
along the outflow axis, and that it has a dependence on the cylindrical radius r which is 
given by a generic Taylor series expansion of the form vj = a + br + cr'^ + . . .. One can argue 
that the flrst order term should be nill, since it leads to an unphysical on-axis "cusp" in the 
distribution. Therefore, the lowest order function that we can consider is quadratic. We will 
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then assume that the radial velocity cross section of the jet is given by 

vj{r) = ^0 (l - ^) , (9) 

where Vq is the central velocity and h the outer radius of the jet beam. It is assumed that 
both Vq and h vary along the jet beam only in scales much larger than the jet width. 

If the jet axis lies at an angle with respect to the plane of the sky (moving away from 
the observer), the profile of an emission line for a point within the jet is given by : 



where 



x'^ + cos^ 



fr = vosm0 1 — . (11) 

L 

In this equation, x is measured across the section of the jet, and z along the line of sight 
(with x, z = 0, on the outflow axis). The line width A^t- includes both the thermal width 
and the turbulent velocity (the instrumental width could be added as well to this term in 
order to reproduce observed results), and is assumed to be independent of position within 
the jet. We also assume that the emitted energy per unit time, volume and solid angle in 
the line % is independent of position (see Canto, Raga & Riera 2003). 

The line intensity at a distance x from the jet axis will then be 

Iy{x) = 2 J,dZ, (12) 



(13) 



where 

COS0 

One can then compute the barycenter velocity 



and the line width 



V,{x)^^vosm<P (l-^) , (14) 



W^^(-) = ^ + |^o^sin^<^(l-^) , (15) 



X 



by carrying out the integrals of equations (4) and (5). 
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Now, from our Fabry-Perot observations we have computed averages of the hne width 
over regions of different characteristic sizes a across the width of the jet. In our model, this 
can be computed as 



<W^> {a, x) 



2a 



W {x)ax = ' ^ 

x—a ^ 



'r—Vo^sm^(l) 
45 



5¥ ~ 3h? 



1-3 



_x 



x'^ a?- 



We can also compute the spatially averaged line center velocity 

<Vc> {a,x) = — I Vc{x')dx' = -fosin 
2 a Jx-a 3 

as well as the square of the standard deviation of the line center velocities 
< Ati^ > (a, x) ^ 



(16) 
(17) 



x 



2 \ 21 



/l2 



2a 



^+'^ 2 16 

[Vc{x')— <Vc> {a,x)] dx' — — i'o"sm' 

x—a 



27 



h 



Equations (16) and (18) define an upper envelope and a lower envelope for the possible values 
that < W'^ > and < Ai>^ > can take as a function of the averaging scale a. These envelopes 
are obtained by setting x — (averaging interval centered on the jet axis), and x — h — a 
(averaging interval touching the outer edge of the jet beam). Setting x = 0, we obtain a 
lower envelope for < W'^ >^/^ and an upper envelope for < Av^ >^/^ : 
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Setting X = h — a, we obtain an upper envelope for < W'^ >^/^ and a lower envelope for 
< Av^ >V2 : 
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The curves for < W"^ >^^'^iow and < W"^ >^^'^up (as well as the ones for < A^;^ >^^'^iow and 
< Av^ >^^^up) meet at the maximum possible value for the averaging length a (i. e., at 
a — h). 

We have drawn the curves given by equations (19-22) on the corresponding plots of 
Figure 8. For each section along the jet we have chosen values of h and Vq sin such that 
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ior a — h (when the upper and lower envelopes join, see above) the model coincides with 

the value of < Av"^ >^/^ corresponding to the largest measured scale length perpendicular to 
the HH 110 axis, which has = a^^max and < Av'^ >i/2= < Aw^ >max^^'^- In other words, 
we choose h = a^^max and fosin0 = (9/4)v^< Av"^ >max^'^- We then choose a value for 
Avt so that for a = h the model prediction coincides with the outer point of the observed 
< W"^ >^/^ vs. a points. The values of /i, Uosin0 and Avt resulting from these fits to the 
different regions along the HH 110 jet are given in Table 2. 

Interestingly, we again find an acceleration along the HH 110 jet, with vq increasing 
from 113 km s~^ (in the A-H region) to 200 km s~^ (in the S region). We can use the 
obtained values of fosin^ to calculate the maximum on-axis value < Vc >max= (2/3)^0 sin 
(obtained by setting x — and a = in equation 17). The resulting values (which are 
also given in Table 2) should be compared with the observations by looking at the on-axis 
values oi < Vc >axis in Figure 7. Prom this figure, we obtain < Vc >axis— —12, —15, —22, 
— 18 and —7 km s~^ for the five chosen regions. If we consider the absohitc vahics of the 
corresponding velocities with respect to the cloud (which has Vhei = +23 km s~^), we obtain 
I < Vc >axis —Vhei\ = 35, 38, 45, 41 and 30 km s^^. For knots A-Q, these values lie within 
20 % of the < Vc >max values obtained from the fits to the < Av'^ >^/^ plots (see figure 8 
and table 2). 

Also interesting is the fact that in order to fit the obcrvcd < W"^ >^/^ (a^.) values, 
we need to include a line broadening of Avx ~ 45 km s^^. If we subtract in quadrature 
the 32 km s~^ instrumental resolution, we obtain an intrinsic broadening Avtm = 35 — 
46 km s^^. These values are substantially larger than the ~ 10 km s~^ expected thermal 
width, so that they therefore probably reflect turbulent motions in the jet flow. 

Looking at Figure 8, we sec that the observed < W"^ >^l'^ and < Av^ >^/^ mostly lie 
on the model predictions for averaging intervals centered on the outflow axis (equations 19 
and 20, respectively). This is not surprising since the wavelets chosen by our knot detection 
algorithm mostly lie close to the HH 110 axis (see figure 3). 

However, for < 5" we see that the observed < W'^ >^/^ and < Av"^ >^^'^ values lie 
substantially above the theoretical curves. In particular, the points corresponding to regions 
that lie further away from the outflow axis (shown with open circles in figure 8) are the ones 
that tend to have larger velocity dispersions. These points can be clearly seen to he in the 
outer, high velocity dispersion envelope which is seen in the cross sections shown in Figure 7 
(also see section 5). We speculate that this envelope corresponds to a turbulent boundary 
layer (which is not present in our simple, analytic jet model) with higher turbulent velocities 
than the central region of the jet beam. 
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7. Conclusions 

We have obtained a Fabry-Pcrot position- velocity Ha cube of the HH 110 flow. Both 
the velocity channel maps and the integrated Ha image show a complex structure of knots 
distributed in a broadening, more or less conical structure. 

Given the description-defying complexity of the structure, we have chosen to carry out 
an anisotropic wavelet analysis of the Ha image, which automatically detects the position 
and characteristic sizes (along and across the jet axis) of the knots. From this analysis, we 
find that the jet first opens out in a cone with a half-opening angle of ~ 10° for distances 
up to 70" from knot A, and has a constant width out to knot R. 

We have then computed the central velocity, the velocity width and the dispersion of 
the central velocity of the line profiles, averaged over ellipses corresponding to the sizes of 
the structures detected by the wavelet analysis. These quantities show that while there is an 
acceleration of the jet out from knot A to knot Q (in the region within 120" from knot A), 
the line width and central velocity dispersion remain relatively constant. A more complex 
behaviour is seen for knots R-S, which could be a reflection of the fact that in this region 
the jet has a curved structure (in the plane of the sky, and therefore probably also along the 
fine of sight) . 

Across the beam of the jet, we see that there is a central, low velocity dispersion region 
with a diameter of ~ 12", surrounded by a higher velocity dispersion envelope with an outer 
diameter of ~ 16". This result could be evidence for the existence of a turbulent boundary 
layer in the outer regions of the HH 110 jet beam. 

Finally, for 5 different regions along the jet beam, we have plotted the mean velocity 
width < >^/^ and central velocity dispersion < At'^ >^/^ as a function of the size scale 
(perpendicular to the jet axis) of the detected structures. We flnd that these plots show 
well organized behaviours. 

Interestingly, we flnd that it is possible to model the observed < >^/^ (ax) and 
< Av^ >^^^ (ttx) proflles with a simple, analytic model of a turbulent jet (or wake). From 
a comparison of the models with the data, we deduce that most of the observed values can 
be explained with a model with an axially peaked mean flow velocity, with superimposed 
turbulent (-hthermal) motions with a velocity dispersion of ~ 40 km s~^. 

However, for characteristic sizes < 5", the observed velocity dispersions lie up to ~ 
25 km s^^ above the model prediction. We flnd that these points correspond to the outer, high 
velocity dispersion envelope of the HH 110 jet (see above). Therefore, in this representation 
we again see the presence of what appears to be a turbulent boundary layer, which is not 
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present in our analytic turbulent jet/wake model. 

Finally, we should note that the acceleration observed along the HH 110 jet does not 
agreee with what one would expect for a turbulent jet. A turbulent jet incorporates low 
momentum material from the surrounding environment, which leads to a progressive decrease 
of the flow velocity as one moves downstream along the jet beam (see, e. g., Raga et al. 1993). 
Therefore, the acceleration observed along the HH 110 jet would favour an interpretation of 
this flow in terms of a turbulent wake, rather than a jet. This is the scenario proposed by 
Raga et al. (1993), who modeled HH 110 in terms of a wake left behind by a jet deflected by 
the surface of a dense cloud, and then pinched off as the incident jet starts to burrow into 
the cloud. 

Of course, another possibility would be to have a turbulent jet ejected with a monoton- 
ically increasing initial velocity. If this increase in the ejection velocity were drastic enough, 
it could overcome the drag due to the entrainment of environmental material, and lead to 
an "acceleration" down the jet axis similar to the one observed in HH 110. 

To conclude, we would like to point out that this work represents a flrst attempt to 
apply wavelet analysis techniques to observations of astrophysical jets. A spectrum of a 
2D image with anysotropic wavelet basis functions is defined in a four- dimensional space 
(with two spatial and two spectral dimensions), and therefore contains a large amount of 
information which is difficult to interpret. The present analysis is based only on a study of 
the spatial and spectral maxima of the wavelet spectrum (i. e., on finding "knots" and their 
characteristic sizes) , and clearly does not use a lot of the information present in the wavelet 
spectrum. 
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Fig. 1. — Ha velocity chanel maps corresponding to various radial velocities obtained from 
the Fabry-Perot data of HH 110. The heliocentric radial velocities are indicated. The 
systemic velocity of the ambient cloud is +23 km s~^. The scale (in arcsec) is given in the 
top left hand plot. The maps are shown with a linear greyscale given (in counts) by the bar. 
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Fig. 2. — The Ha image map (left panel), central radial velocity map (central panel) and the 
line widths (FHWM, right panel) obtained from the position-velocity cube. We have rotated 

the maps by 13° so that the outflow is more or less parallel to the y-axis. The greyscales 
of the map (linear, in counts) and of the central radial velocity and line width (linear, in 
km s~^) are given by the bars on top of each plot. 
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Fig. 3. — Left panel: contour plot of the Ha image, where we have labelled the knots to help 
with their identification. Right panel: the positions and characteristic sizes of the structures 
along the HH 110 jet. The crosses show the positions of the maxima Xk of the wavelet 
transform obtained for the different values y along the jet. The characteristic sizes a^^k and 
of these maxima are shown as error bars centered on the positions of the maxima. The 
distances x (across the jet axis) and y (along the jet axis) are measured from the peak of 
knot A. We have rotated the maps by 13° so that the outflow is more or less parallel to the 
y-axis. 
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y ["] 

Fig. 4. — The < > arc plotted as a function of position y along the HH 110 jet. A straight 
line that forms an upper envelope to the < ax > values (in the y < 70" region) has been 
drawn (see the text). 
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Fig. 5. — The Qy^k are plotted as a function of ax,k for different regions along the HH 110 
flow (see Table 2). For regions 0-Q and R1-R2 the results from linear fits to the points are 
drawn as solid lines. 




Fig. 6. — The average line centers (top panel), line widths (central panel) and standard 
deviations (bottom panel) as a function of position y along the HH 110 jet. All points with 
different positions x across the jet and with different and ay are plotted. 
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Fig. 7. — The line centers (left), line widths (center) and standard deviations (right) are 
plotted as a function of position x across the jet beam for the 5 different chosen regions 
along the HH 110 flow, xq is the central position of the jet cross section (which varies along 
the jet, as the beam has side-to-side excursions, see figure 3). 
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Fig. 8. — The average standard deviations of the hne center (left panels) and the hne widths 
(right panels) obtained from the data cube are plotted as a function of the size scale < ax > 
across the jet flow for the 5 chosen regions along the jet. The triangles correspond to points 
which lie in the central region of the flow (with |a;| < 3" for the A-H region, |a;| < 7" for the 
1-N region, < 5'.'5 for the 0-Q region, and \x\ < 9" for the R1-R2 region), and the open 
dots correspond to the outer region of the jet beam (see section 6). The curves correspond to 
the predictions of the analytical model, which have been fitted independently for each region 
along the HH 110 jet (see section 6). In the left panels, the curves correspond to < Av'^ > 



(solid line) and < Av'^ >\!Z (dashed line). In the right panels, the curves correspond to 
< >\!p (sohd line) and < >J/^ (dashed hne). 
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Table 1. Journal of Fabry-Perot observations 



Observations 


Telescope 


ESO 3.6m 




Equipment 


CIGALE @ Cassegrain 




Date 


1997, January, 4-6 




Seeing 


0'.'4 O'.'S 


Calibration 


Neon Comparison light 


A 6598.95 A 


Fabry-Perot 


Interference Order 


793 @ 6562.78 A 




Free Spectral Range at Ha 


380 km s-i 




Finesse^t Ha 


11 




Spectral resolution at Ha 


9400*' 


Sampling 


Number of Scanning Steps 


24 




Sampling Step 


0.35 A(16 km s"!) 




Total Field 


170" X 170" (256 X 256 px^) 




Pixel Size 


0.91" 


Detector 




IPCS (Thomsom tube) 



''Mean Finesse through the field of view 

''For a signal to noise ratio of 5 at the sample step 



Table 2. Model fits to the HH 110 jet 



Region Dmin, Umax h sin VO <Vc>max Avt ^VT,int 

["] ["] [km s-1] 



A-H -11, 53 15 65 113 43 52 46 

I-N 53, 85 18 65 113 43 51 45 

O-Q 85, 121 14 83 145 55 48 41 

R1-R2 121, 136 19 95 166 63 45 37 

S 144, 171 16 115 200 77 43 35 



Note. — Vmin, Umax give the position of the chosen regions along the HH 110 jet, 



where y = corresponds to the peak of knot A. h is the outer radius of the jet beam 
and VO sin </> is the central radial velocity of the model fit to the data (see section 6) . 

The deprojected velocity vq was computed assuming an inclination angle = 35° 
with respect to the plane of the sky. < Vc >max is the maximum on axis value for 
the baryeeiiter of the line profile computed from the model fit. Avt is the turbu- 
lent+thermal+instrumental line width necessary to fit the observed < >i/2 vs. a 
(see section 6). AvT,int is the intrinsic broadening, computed substracting in quadra- 
ture the instrumental width from Avt , where Avi^s = {FHWMi„s)/{2{ln2)'^/'^). For 
an instrumental FWHM = 34 km s~^, a Avi„s = 25 km is obtained. 



